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Spectral properties of SrVOs 



O 
(N 



\0' 



Himpyo Lee,-'^ Kateryna Foyevtsova/ Johannes Ferber/ Markus Aichhorn,^ Harald O. Jeschke/ and Roser Valentf-'^ 

^Institut fiir Theoretische Physik, Goethe- Universitdt Frankfurt, 
Max-von-Laue-Strafle 1, 60438 Frankfurt am Main, Germany 
'^Institute of Theoretical and Gomputational Physics, TU Graz, Petersgasse 16, Graz, Austria 

(Dated: March 27, 2012) 

We present a combination of local density approximation (LDA) with the dynamical cluster 
approximation (LDA+DCA) in the framework of the full-potential linear augmented plane-wave 
method, and compare our LDA-I-DCA results for SrVOa to LDA with the dynamical mean field 
theory (LDA-I-DMFT) calculations as well as experimental observations on SrVOs. We find a 
qualitative agreement of the momentum resolved spectral function with angle-resolved photoemission 
spectra (ARPES) and former LDA+DMFT results. As a correction to LDA+DMFT, we observe 
more pronounced coherent peaks below the Fermi level, as indicated by ARPES experiments. In 
addition, we resolve the spectral functions in the Kq — (0, 0, 0) and Ki = (tt, tt, tt) sectors of DCA, 
where band insulating and metallic phases coexist. Our approach can be applied to correlated 
compounds where not only local quantum fiuctuations but also spatial fluctuations are important. 
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I. INTRODUCTION 

The development of reliable numerical tools for the de- 
scription of the electronic structure of correlated com- 
pounds is one of the most challenging tasks in the con- 
densed matter community. As an example, transition- 
metal perovskites with partially filled t2g orbitals are pre- 
dicted to be conventional metals in the framework of one- 
electron approaches like density functional theory (DFT) 
in the local density approximation (LDA). Nevertheless, 
a few perovskite families show a markedly different be- 
havior; SrVOs and CaVOa are correlated metals with sig- 
nificant mass enhancement and LaTiOs displays features 
of a Mott insulator In all compounds, this anoma- 
lous behaviour may be caused by correlations resulting 
from Coulomb repulsion effects. Therefore, progress on 
methods including correlation effects beyond DFT is very 
desirable. 

Dynamical mean field theory (DMFT)2r— takes local 
quantum fluctuations fully into account but the momen- 
tum dependence of the self-energy is neglected. This 
method has been developed over the last twenty years 
and successfully describes the metal to Mott insulator 
transition in frustrated systems^^^'' and non Fermi-liquid 
behavior in multi-orbital systemsiir— to mention a few 
examples. On the other hand, it cannot describe such 
phases as spin density wave or d-wave superconductiv- 
ity due to its lack of spatial correlations. In order to 
overcome these problems the first extension to single- 
site DMFT are multi-site approaches in which the short- 
range correlations are exactly considered within a clus- 
ter—^. Implementations of these approaches are the 
dynamical cluster approximation (DCA)i^ii^ or the cel- 
lular DMFT~. These approaches capture the spin den- 
sity wave formation indicated as a band insulatorii'i^ as 
well as Mott transitionsi^Ti^i. Recently, other implemen- 
tations including long-range correlations have been de- 



veloped by considering a perturbation expansion where 
nonlocal contributions are obtained from the two-particle 
vertex functions^"—. 



Recent progress towards a realistic description of corre- 
lated systems is the combination of LDA with DMFT— 
(LDA-I-DMFT). While this approach has proven to be 
quite successful for the description of spectral properties 
of transition metal oxides^i^"— and the newly discov- 
ered iron-based superconductors^—, effects originating 
from spatial fiuctuations remain inconclusive. Attempts 
to include short-range spatial fiuctuations have been done 
in the context of the spin-Peierls system TiOCl -where 
pairing correlations are important- within an N*'' or- 
der muffin tin orbital (NMTO) approach combined with 
DCA^^ as well as NMTO combined with a variational 
cluster approach (VCA)i^. In this work, we present an 
alternative approach where we extend a newly developed 
implementation of the LDA+DMFT approach^^ in the 
context of the full-potential linearized augmented plane 
wave (FLAPW)^ method by including spatial fluctua- 
tions within DCA (LDA-I-DCA), and we investigate the 
spectral properties of SrVOa as a test case. 



The paper is organized as follows: in section |TT1 we 
describe our LDA-I-DCA implementation with a weak- 
coupling continuous-time quantum Monte Carlo (CT- 
QMC) algorithm""^"^^ for multi-orbital systems with mul- 
tiple sites. In section Hill we present results for SrVOa 
within LDA-I-DCA with a cluster of two sites and com- 
pare them with single-site LDA+DMFT calculations as 
well as experimental observations and in section IIVI we 
summarize our findings. 
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II. THEORETICAL FRAMEWORK 
A. LDA+DCA in the APW framework 

In this work, we extend a recent implementation of 
LDA+DMFT^i to LDA+DCA which includes short- 
range spatial correlations. We first shortly review the 
projection operators within the WIEN2K code^. The 
local atomic-like Wannier orbital functions inside an ap- 
propriate energy window W can be expanded over the 
Bloch basis set as 



where E^™?(K, ia;„) is calculated by the continuous 
time Quantum Monte Carlo (CT-QMC) cluster solver 
and the Dyson's equation. Calculating the DC correction 
is not possible exactly, but some approximate expressions 
have been introduced. Here, we use as double counting 
correction 



(7) 



(1) 



where a indicates the correlated atom, v is the band in- 
dex, a the spin index, and m the orbital index. Here, 
hAk i^) is the Bloch eigenfunction in the augmented plane 
wave basis and the correlated orbital \x"f) is given 
as \Xm'^} = Wf''^ {Ei)Yjn) within the muffin tin sphere, 
where Ei are chosen linearization energies, '"^ is the 
radial wave function, and Yj"^ is the spherical harmonic 
function. The orthonormalized projector operators for 
the DMFT and DCA self-consistent equations are calcu- 
lated by 

(2) 



where U' = U — 2J , U is the onsite Coulomb interac- 
tion, J is the Hund's coupling and Nc and denote the 
number of total occupied states and spin-resolved occu- 
pied states in the correlated orbitals, respectively^. The 
local cluster Green's functions are given as 

= J2 C:^(K + k)G-,,(K + k,ic.„)P,"/;;'(K + k), 



(8) 



where the summation over k is calculated in each Bril- 
louin zone sector. The LDA-I-DCA self-consistency con- 
dition states that this local cluster Green's functions, Eq. 
([H]) , have to be equal to the impurity Green's functions as 
calculated by CT-QMC. The DMFT update of the Weiss 
field is given by the Dyson's equation as 



where 0(k, a)Xn' is the overlap function which is given [^'^o^^ (k, ' ^ E:;X?(K, za;„)+ [G"J°^, (K, zc.„)] 



as 



0(k,a)^i, = E (x'^kli'^JiKM^n- (3) 



For the LDA-I-DCA self-consistency procedure, the lat- 
tice Green's function is given as 



G^^,,(K-Fk,zw„) 



1 



K+k,i 



S]-,,(K + k, 

(4) 

where we have defined k= K + k with K being the cluster 
momenta and k running over each Brillouin zone (BZ) 
sector. cOn is the Matsubara frequency, is the chemical 
potential, e^^j^ ^ are the Kohn-Sham (KS) eigenvalues, 

and EJj,/(K -I- k, is the lattice self-energy which is 
calculated as an expansion of the cluster self-energy over 
the Bloch basis set: 

E^_,,(K + k,iL^„) 

P,";^' (K + k) AE^X?(K, iu^n)Py,, (K + k). 



(9) 



B. Many-body interactions and CT-QMC 
algorithm 

In order to describe the electronic behavior of SrVOs 
one has to consider the multiorbital Hubbard Hamilto- 
nian where the interaction term is given by: 



E 



a.m.m 



(5) 

From the self-energy we need to subtract the contribu- 
tion to correlations that is already included in the LDA 
calculation, commonly called double counting (DC) cor- 
rection. 



AI]--?(K,^c.„) 



y^CT, imp 

771, rn' 



(K,iw„) - E 



dc 



(6) 



+ E 

m m<Ti,cr 

and m,n denote the t2g orbitals. In order to solve this 
model we employ a weak-coupling CT-QMC algorithm. 
While the weak-coupling CT-QMC algorithm can eas- 
ily treat a multiple number of sites in the cluster, it is 
difhcult to deal with the full rotationally invariant form 
of the interaction Hamiltonian due to the fermionic sign 
problem, in contrast to the strong-coupling CT-QMC al- 
gorith m^^'^" . Therefore, in what follows we shall consider 
a simplified Hubbard model where the spin-flip and pair- 
hopping terms in Eq. pH]) are neglected (J' — 0). 

The main idea of the weak-coupling CT-QMC method 
is to divide the total action S into an unperturbed term 
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So and the interaction term / which is expanded in a 
Taylor series. The partition function is rewritten as 

Z = Y.Z,t^ j dT,.--drJv[cA 

(ni,^(Ti)ni.;(ri) • • • n,,^(Tfe)n,- ;(Tfc)), 

where (n-i,t(Ti)rii;4(ri) • • ■ rii^-l.(rfc)n;;^4,(Tfc)} is deter- 
mined by the non-interacting Green's function and 
Wick's theorem, fc is the perturbation order, Zq — 
Tr(re~'^°) corresponds to the unperturbed term and /, 
and Tfc are randomly sampled. / is given as 

I ^UPNM{2M -1), (12) 

where (3 is the inverse temperature, N and M are the 
number of sites and the number of orbitals in the cluster, 
respectively, and U is one of U, U' or U' — J depending 
on the operators considered in the random walk in the 
average (...) in Eq. (|lip . The impurity Green's functions 
are calculated by numerically averaging Eq. pT|) . 

III. RESULTS 

SrVOs, which is thought to be a prototypical para- 
magnetic correlated metal with intermediate electron- 
electron interactions, has served in the past as a test- 
ing ground for numerous newly developed LDA-f DMFT 
approaches'*^"— In SrVOa, the V 3d-orbitals are 
split by the crystal field into triply degenerate t2g and 
doubly degenerate eg states. The LDA calculations show 
that the degenerate t2g states of V form bands cross- 
ing the Fermi level which are well separated from the Cg 
bands. 

For our calculations on SrVOa, we chose the energy 
window W from -1.35 eV to 2.0 eV for the t2g or- 
bitals which can then be effectively described by the 
degenerate three-orbital Hubbard model in Eq. ([TU|. 
We first reproduced the results of LDA+DMFT from 
Ref. HH, considering a temperature T = 0.1 eV and the 
same Coulomb interaction U = 4.0 eV, and Hund's rule 
coupling J = 0.65 eV. In a next step, we extend the 
LDA-I-DMFT solution to LDA+DCA with two sites in 
the cluster iV = 2. Within the DC A method, = 2 
implies that we have two BZ sectors and the self-energies 
in the BZ sectors are constant. 

The BZ of SrVOs has cubic symmetry and the 
self-energies in the cluster momenta Ko = (0,0,0) 
and Ki = (7r,7r,7r) are calculated in the BZ sec- 
tors shown in Figs. [Ha) and (b). In real space, 
the on-site and nearest neighbor-site Green's func- 
tions are G_r=o(«w„) = i(G'Ko(«w„) + Gki(«w„)) and 
G'fl=i(iw„) = i(G'Ko(«w„) - GKi(ia;„)), respectively. 
Here, the DCA formalism with = 2 for cubic lattice 
has been clearly presented in Ref. [5l|. Both on-site and 
nearest neighbor-site Green's functions are inserted into 
the CT-QMC impurity solver, and the LDA-I-DCA self- 
consistency is satisfied by Eqs. (|S]) and 




FIG. 1. (Color online) One quadrant of the Brillouin zone 
of SrVOs. (a) and (b) represent the Ko = (0, 0, 0) and Ki = 
(tt, tt, tt) Brillouin zone sectors, respectively. Other quadrants 
follow from symmetry. 
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FIG. 2. (Color online) (a) The density of states p{Ld) of 
SrVOa calculated within the LDA-I-DMFT and LDA-fDCA 
approaches, with U = 4.0 eV, J = 0.65 eV and T = 0.1 eV. 
The density of states p{uj) for LDA-I-DCA is calculated by 
p{uj) = i(A(Ko,w) + A{Ki,Lj)), with Ko = (0,0,0) and 
Ki — {n,n,n) sectors, (b) and (c) The spectral functions 
A(K, cj) obtained from LDA and LDA-I-DCA for the Ko and 
Ki sectors by Eq. ([8}, respectively. All the density of states 
and spectral functions are normalized to one. 
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Fig. [2ja) shows the density of states p{uj) of the vana- 
dium orbitals obtained within the LDA+DMFT, 
p{uj)= A{uj) and LDA+DCA, p{uj) = ^{A{Kq,u) + 
A(Ki,w)). Here, the spectral function A(K,iuj) is given 
as 



DMFT 



1_ 



A{K,uj) = ImGK(w), 



(13) 



and an analytical continuation of the impurity Green's 
functions is performed through a maximum entropy 
method. Our LDA+DMFT resuhs obtained by both 
the weak-coupling CT-QMC^i^ as well as the strong- 
coupling CT-QMC algorithms from the ALPS code^ 
agree with former LDA-I-DMFT calculations^i^i^^"^. In 
Figs. IHb) and (c) we present the spectral functions for 
the Kg — (0,0,0) and Ki = (tt, 7r,7r) sectors within 
LDA and LDA-I-DCA. The new features obtained in 
LDA-I-DCA are a broad peak around 1.5 eV and a co- 
herent peak around 0.2 eV below Ep. LDA results (see 
Fig. Mjo)) as well as most former LDA-fDMFT results 
(see also Fig. [2] (a)) don't exhibit neither a broad peak 
nor a clear coherent peak below Ep. Recent angle re- 
solved photoemission (ARPES) experimental, have ob- 
served, in fact, a broad peak around 1.5 eV and a co- 
herent peak around 0.4 eV below Ep (Fig. 3(b) in 
Ref. |53I ). We suggest that the better agreement of 
the LDA-I-DCA with ARPES observations is a conse- 
quence of the inclusion of short-range spatial correlations. 
Fig. HIc) shows that the coherent and broad peaks below 
the Fermi level are caused by the distinct spectral weights 
in the Kq — (0,0,0) and Ki = (tt, 7r,7r) sectors. These 
two sectors also show respectively metallic and band in- 
sulating behavior reminiscent of the LDA results in these 
sectors (FigHJb) and Ref. HI. 

In Figs. [3] (a) and |3] (b), we show the momentum re- 
solved spectral functions calculated from Eq. Q and 
Eq. (jSl). The analytical continuation of the self-energy 
S(K,jcj„) is performed by the maximum entropy ap- 
proach with subtraction of the Hartree-Fock term^S. In 
view of the ill-posed problem of the analytical contin- 
uation of the self-energy^, we also compared the DOS 
obtained from integration of the spectral functions with 
those in Fig. [5] (a) and found a reasonable agreement. 
We also compare LDA-hDMFT to LDA-I-DCA resuhs. 
One can observe some redistribution of momentum re- 
solved spectral weight between the LDA+DMFT and 
LDA+DCA resuhs. In Figs. [3] (c) and O (d) we plot 
the LDA+DMFT and LDA+DCA spectral functions re- 
spectively, in the region between ky = —0.6-^ and 0.4^ 
at kx ~ 0.0 and kz = 0.32-^ (a is the lattice con- 
stant) in order to directly compare our calculations to 
the ARPES resuhs (Fig. 1 (a) of Ref^). In agree- 
ment with ARPES experiments, both LDA+DMFT and 
LDA+DCA show dispersive features around -0.7 though 
they are more pronounced in the LDA+DCA calcula- 
tions. Also, the LDA+DCA calculations reproduce the 
small peak observed around -0.2eV. These results account 
for the renormalization of the bands due to electronic cor- 
relations. 
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FIG. 3. (Color online) (a) and (b) Spectral functions obtained 
within (a) LDA+DMFT and (b) LDA+DCA for the vana- 
dium t2g bands, (c) and (d) Spectral functions in the regime 
between ky = -0.6^ and 0.4f at fc^ = 0.0 and = 0.32^ 
within (c) LDA+DMFT and^) LDA+DCA for the vana- 
dium t2g bands. 



Finally, our estimation of the mass enhancement is 
m*/m « 1.7 ± 0.3 within LDA+DMFT and 1.6 ± 0.5 
within LDA+DCA at Kq = (0,0,0). These values are 
obtained from 



m* /m « 1 



9ImS(iLLin) 



>o+ 



(14) 



where the derivative is extracted by fitting a third- 
order polynomial to the lowest four Matsubara frequen- 
cies^. Note that our LDA+DCA estimates give a slightly 
smaller mass enhancement than LDA+DMFT estimates 
with a larger error bar. Both sets of results are in accor- 
dance with ARPES estimates of m*/m w 1.8 ± 0.2^. 



IV. CONCLUSIONS 

In conclusion, we have presented an implementa- 
tion of the LDA+DCA method within the linear aug- 
mented plane-wave framework. We have compared 
our benchmark results on SrVOs, which is modeled in 
terms of a three-band Hubbard Hamiltonian, with ear- 
lier LDA+DMFT calculations as well as experimental 
data. Since the LDA+DCA approach considers both lo- 
cal quantum as well as short-range spatial fluctuations, it 
offers a more complete description of correlated materi- 
als compared to the LDA+DMFT approach, where only 
local quantum fluctuations are taken into account. 

Unlike the LDA+DMFT, the LDA+DCA approach re- 
produces both coherent and broad peaks for SrVOa be- 
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low the Fermi level, observed in angle integrated photoe- 
mission experiments. The analysis of the spectral func- 
tions at Kq = (0,0,0) and Ki — (tt, 7r,7r) reveals the 
source of these peaks. While the broad peak is due to 
the spectral function in the Kg = (0, 0, 0) sector, the 
coherent peak has its origin in the spectral function at 
the Ki = (tt, tt, tt) sector. We also observe a metallic 
and a band insulating state at the Kq = (0,0,0) and 
Ki = (tt, tt, tt) sectors, also present in the LDA results. 

In summary, we believe that the presented LDA+DCA 
approach is very promising and can be applied to a large 
variety of multiorbital correlated compounds at different 



fillings. 
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